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Abstract 

We expound in detail a method frequently used to reduce 
the Dirac equation in D-dimensional {D > 4) spherically 
symmetric spacetimes to a pair of coupled partial differ- 
ential equations in two variables. As a simple application 
of these results we exactly calculate the quasinormal fre- 
quencies of the uncharged Dirac field propagating in the 
I?-dimensional Nariai spacetime. 



1 Introduction 

Recently in many research lines of theoretical physics the 
models in which the spacetime has more dimensions than 
the four dimensions observable in our daily experience have 
been studied extensively. The most analyzed models are 
those related to string theory [T] . Also the scrutiny of the 
properties and solutions of higher dimensional general rel- 
ativity has attracted a lot of attention (see Ref. [2] and ref- 
erences therein) . In several of these research lines we need 
to know the classical properties of the higher dimension- 
al spacetimes to examine different phenomena. Therefore 
the investigation of these classical properties is an active 
research field. 

To analyze the classical properties of a given spacetime 
a common method is to use a field as probe [3], [4]. Thus 
in the past several scattering phenomena of classical fields 
were studied, in order to know how to calculate the physi- 
cal parameters of the spacetime from the measured values 
of the physical quantities corresponding to the classical 
field. 

The quasinormal modes (QNMs) are solutions to the 
equations of motion for a classical field that satisfy the ra- 
diation boundary conditions that are natural in the space- 
time in which the field is propagating [3], [3]. For example, 
in asymptotically flat black holes the boundary conditions 
of the QNMs are that the field is purely ingoing near the 
event horizon and purely outgoing near infinity [3]. For 
asymptotically anti-de Sitter black holes we impose the 
boundary condition that the field vanishes at infinity and 
is ingoing near the event horizon. 
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It has been shown that the QNMs are a useful tool to 
calculate the physical parameters of a spacetime [3], 
Hence if we know the quasinormal frequencies (QNF) of 
a classical field we can infer the values of several physi- 
cal quantities of the spacetime such as its mass, charge, 
and angular momentum [Sj. Furthermore it has been pro- 
posed that the QNMs encode some information about the 
quantum properties of the black holes [S] . 

To compute the QNF of a classical field in a given space- 
time the usual procedure is to reduce the equations of mo- 
tion for the field to a radial ordinary differential equation 
(assuming a given dependence on the angular variables and 
a harmonic time dependence) and impose to the radial 
function the boundary conditions of the QNMs. 

Also notice that the reduced form of the equations of 
motion is useful (and sometimes necessary) to study many 
other classical or semiclassical phenomena. Thus we be- 
lieve that at present time the understanding of the sepa- 
rability properties of the equations of motion for classical 
fields in higher dimensional curved spacetimes must be a 
relevant part in the education of a physicist. 

Motivated by these theories that assume a number of 
spacetime dimensions greater than four, the separability 
properties of the equations of motion for several classical 
fields were studied in higher dimensional backgrounds. It 
was found that many of the well known results that are 
true for four-dimensional spherically symmetric spacetimes 
extend to D-dimensional {D > 4) spherically symmetric 
spacetimes. 

For example, the reduction of the equations of motion for 
Klein-Gordon, electromagnetic, and gravitational pertur- 
bations to ordinary differential equations which is true in 
four-dimensional uncharged spherically symmetric space- 
times j4 , also is valid in Z3-dimensional uncharged spher- 
ically symmetric spacetimes, as showed in Refs. [6], [7]- 
Moreover for the coupled gravitational and electromag- 
netic perturbations the reduction of the equations of mo- 
tion to Schrodinger type equations which is true in four-di- 
mensional charged spherically symmetric backgrounds also 
is valid for _D-dimensional charged spherically symmetric 
backgrounds [8]. 

Although the study of the classical dynamics of fields in 
curved spacetimes (in four and D dimensions) is focused 
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on boson fields [3], [4], [6], [7], [E] mainly on gravitational 
perturbations, we believe that the understanding of the 
classical dynamics of the fermion field in D-dimensional 
spacetimes is of great value, because the Dirac field some- 
times behaves in a different way that the boson fields. For 
example, it is a well known fact that in the four-dimen- 
sional Kerr black hole the fermion field does not show su- 
perradiant scattering [9] , [10] , unlike to boson fields [4] . 

For the Dirac equation, its separability properties in D- 
dimensional spherically symmetric spacetimes were previ- 
ously studied in Refs. [11], [12]. In these papers is shown 
that the Dirac equation reduces to a pair of coupled partial 
differential equations in two variables 

Owing to the past and future applications of the reduced 
system of partial differential equations for the Dirac equa- 
tion in Z)-dimensional spherically symmetric spacetimes, 
we believe that the method used in Refs. [11] to reduce 
the Dirac equation to a pair of coupled partial differential 
equations deserves a detailed exposition, because this ac- 
count may be practical and useful. Here we present the 
method in more detail than in the original references, that 
is, in the present work we explicitly write some mathe- 
matical steps omitted in Refs. [11], [12] (see also [T^, [16], 

Notice that Section [2] is not an exhaustive review of the 
previous work on the dynamics of fermion fields in spheri- 
cally symmetric spacetimes. Also observe that in this pa- 
per we do not consider in detail the mathematical prop- 
erties of I?-dimensional spinors, these can be studied in 
many books and papers (see for example Refs. [18], [19] . 
[20]!. We only write the essential properties of the spinors 
necessary to make the reduction of the Dirac equation to 
two coupled partial differential equations that we shall ex- 
pound in Section 121 

Recently the exact computation of the QNF for several 
higher and lower dimensional spacetimes has attracted a 
lot of attention, see [lT]-[3ni for some references in which 
an exact calculation of the QNF was carried out. As many 
exactly solvable models in theoretical physics, we believe 
that these examples are useful models and it is possible 
that they play a relevant role in future research. 

The Nariai spacetime is a vacuum solution to the Ein- 
stein equations with positive cosmological constant [40] , 
This spacetime is a simple solution to the field equations 
of general relativity. Owing to this simplicity of the Nar- 
iai solution, it is possible to calculate the values of several 
physical quantities in exact form. For example, for this 
spacetime in Refs. [39. were computed exactly the values 
of the QNF for Klein-Gordon fields and tensor type grav- 
itational perturbations. 

Furthermore in the £)-dimensional charged Nariai space- 
time [8^ , [40j , the QNF for coupled gravitational and elec- 
tromagnetic perturbations were calculated exactly in Ref. 
[38] . To our knowledge the result of the previous reference 

^Notice that for Dirac field some results valid in four-dimension- 
al rotating black holes have been extended to rotating black holes 
in higher dimensions, see Refs. |13l for an incomplete list of refer- 
ences. For a review of the recent work on the separability properties 
for the equations of motion for several fields in higher dimensional 
spacetimes see Ref. |14l . 



in the charged Nariai spacetime is the only exact calcula- 
tion of QNF for coupled electromagnetic and gravitational 
perturbations in higher dimensions. 

As an application for the reduced system of differential 
equations obtained in Sectionj^Jfor Dirac field moving in D- 
dimensional spherically symmetric spacetimes, we exactly 
calculate the QNF of this field in _D-dimensional Nariai 
spacetime [40]. These values of the QNF for Dirac field 
extend those already published in Refs. [38], [39] . 

In this paper we assume that the reader has a working 
knowledge of general relativity and differential geometry. 
Furthermore, in the following sections we use Einstein's 
sum convention and understand sum on repeated indices 
(Latin and Greek indices), unless we explicitly state that 
in a given formula we do not understand sum on repeated 
indices. 

The paper is organized as follows. In Section [2] we 
present in detail the method of Refs. [11] (see also [15], 
[l6] . [17] ) that reduces the Dirac equation in D-dimen- 
sional spherically symmetric spacetimes to a pair of cou- 
pled partial differential equations in two variables. Using 
these results in Section [3] we exactly calculate the QNF 
of the Dirac field propagating in the D-dimensional Nariai 
spacetime. Finally in Section |3] we discuss some related 
facts. 

2 Dirac's equation in D-dimen- 
sional spherically symmetric 
spacetimes 

As is well known in D-dimensional spherically symmetric 
backgrounds the Dirac equation 

iJJip = rnip, (1) 

reduces to a pair of coupled partial differential equations 
in two variables [TT] , [T^] , [H] , [H] , [T7] . In this section we 
describe in detail the method of Refs. [TT] frequently used 
to get this result. For a different method see Refs. [T^ . 

Here we shall consider two D-dimensional spacetimes M 
and M whose metrics g^j/ and g^, are conformal, that ig 

where £7 is a function of the coordinates. We point out 
that if the symbols p^ip, m, and p, ip, ifi denote the Dirac 
operator, the Dirac spinor, and the mass of the Dirac field 
corresponding to the spacetimes M and M, respectively, 
then the following relations are satisfied [11] , [15] , [16] , [17] 

fiP = (3) 
m = Qm. 

It is well known that the previous results can be gener- 
alized when there are gauge fields [llj . but in this paper 

^Notice that Greek indices stand for coordinate indices, whereas 
the Latin indices stand for frame indices. 
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we do not analyze this extension. As in formulae ([3]), in 
the rest of the present section a tilde stands for quantities 
corresponding to the spacetime with metric g^j,^. 

One way to obtain the results of formulae ([3]) is the 
following. When the metrics g^v and g^^ are conformally 
related as in formula ([2]) and we define the basis of one- 
forms e° such that = e^ej^ryab and the basis e° such 
that = el^e^rjab, where rjab = v"'' = diag(l, -1, . . . , -1) 
is the Minkowski metric [41], [42], we find that the one- 
forms e" and e" satisfy 



(4) 



To get the relation between the connection one-forms 
a;"j and ti^j, corresponding to the basis of one-forms and 
e", respectively, we recall that the connection one- forms 
uj°'fj are determined by the first Cartan structure equation 



de"^ 



(5) 



where the symbol A stands for the wedge product [H] , [H] . 

Using Statemenl[f| 6.1.6.1 of Ref. [41], we see that the 
one-forms tD^j, and lj^j, are related by 



u}ab = ^^ab + -^ieb{^)ea- ea{fl)eb), (6) 

where ea{^) denotes the action of the vector = rjabs'' 
on the scalar function fl. From expression ^ we get 

i^abc= '^^^ + -^{eb{^)riac- ea{^)ribc), (7) 

where Uabc = ^abiec) and similarly for uiabc- 

Thus if the symbol Vc stands for covariant derivative of 
a spinor, that is, [18], [42] 



(8) 



where stands for Z3-dimensional gamma matrices that 
satisfy [IS], [H, [ID] 



+ 77" = 27/'^^ 



(9) 



and we also observe that 7"^ = 7" . Taking into account the 
previous results and if we take ip = ri'^ip then it is possible 
to show that the Z?-dimensional Dirac operator 



(10) 



•^Statement 6.1.6.1 of Ref. |41j : If any set of 2-forms Ac is given 
and {0"^} is a dual frame basis then there exists a unique set of 1- 
fomis Xab such that 

Aa + '^Xab^e'' =0, Xab = -Xba- 

b 

The 1-forms Xab can be expressed by the formula 

Xab ^ ' \abc0 : )Cabc — ^{-^abc ^bac ^cab)) 

c 

where A^ab ^^re the coefficients in the decomposition 

= — ^ ^ ^cab^"^ ^ ^cab = ~^cba- 



transforms into [TT], [15], [16], [17] 



flCc + - {i^abc^ + ea{^)ribc 



Sb{n)r,ac)ri']^'i' (11) 



Taking s = (D — l)/2 in the previous formula we finally 
get the result 

= (12) 

This expression and -0 — fi*-'^^^-'/^?/' are the first two re- 
sults of formulae The result given in expressions ^ 
for mass m immediately follows from the previous formulae 
for p-ip and ip. 

In the following paragraphs we study the Dirac equa- 
tion in the D-dimensional spherically symmetric space- 
times with coordinates {t, r, where j = 1, 2, . . . , {D—2), 
and whose line elements we write in the form 



(13) 



where F{r), G{r), and H{r) are functions only of the co- 
ordinate r, the symbol dS|,_2 stands for line clement of 
a (D — 2)-dimensional invariant base space which depends 
only on the coordinates 4'i. 

To simplify the Dirac equation in a spacetime whose 
metric takes the form we factor out the function 

H{ry in the line element and define 



wher^ 



dS^ 



ds'^ 



(14) 



-g^Ax^dx'' = -j^dt' - —dr' - dSl5_2- (15) 



Next we use the results (jSj) to find the relation among 
the quantities ■0, pil^i ™j and ip, pi^, rh, corresponding 
to the spacetimes with line elements ds^ of formula 
and ds^ of expression (fT5]) . We point out that in the line 
element ([T5]) the first two terms depend only on the two 
coordinates t, r, and the term dS]|,_2 depends only on the 
{D — 2) coordinates 4'i. 

As a basis of one-forms for the spacetime with metric 
(Ull) we choose [11], [15], [16], [17] 

g* = g* (r, t) = /(t) (r, t)dt + /(,) (r, t)dr, 

= i-{r, t) = (r, i)di + g^r){r, t)dr, (16) 

? = e^(</'fe) = /i;('^feM<^'' 



*We shall write the functions F{r), G{r), and H{r) simply as F, 
G, and H, respectively. In general, we shall use a similar convention 
for the functions that we shall define in the following paragraphs. 
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where i,j,k— 1,2, — 2. We prefer this basis because 
many of the connection one-forms are equal to zero, for 
example, the connection one-forms ojtj and Qrj- 

It is convenient to observe that in D even dimensions the 
gamma matrices are square matrices of dimension 2^/^ x 
2^/2 whereas in D odd dimensions these are of dimension 
2(r>-i)/2x2(i5-i)/2 [ig^ ^ ^ Thus if in the ZJ-dimen- 

sional spacetime with metric (fT5|) we use the representation 
of the gamma matrices [11] , [15] , [16] , [17] , [19] , 



option, also the quotients of hyperbolic spaces are possible 

[Ml, [^, [13j. 

Taking the spinor ijj of the spacetime with line element 
(fTS)) in the form 



(21) 



(Jl ® l2(D-2)/2, 
= i(T2 ® l2(D-2)/2 , 

= (73® 71, (17) 

= 0-3 «)72, 
= (73 (g) 73, 



7t = tJl (g) (To ® tTo (X) (To ■ • • 

7r = ?(T2 ® (To (8) (To ® (To • • 

71 = *(T3 ® (Tl (8) (To ® (To • • 

72 = *(T3 ® (T2 (8) (To ® (To ■ ■ 

73 = «(T3 (T3 (Tl (g) (To • • 

7D-2 = • • • = (Ta (g) 7D-2, 

where I2(-d-2)/2 is the identity matrix of dimension 
2(-D-2)/2 ^ 2(-D-2)/2^ the symbol (g stands for direct prod- 
uct [H], [H], 7i, 72, 7D-2 are a representation of 
the gamma matrices for a (Z? — 2)-dimensional space with 
signature (—.....—) and 



(To = 



where ^2_d {f, t) is a two-spinor on the spacetime with line 
element ds2D of formula (|20|) and the functions x{^i) sat- 
isfy 

Vd^X = ^X, (22) 
that is, X and k denote the eigenfunctions and eigenvalues 
of the Dirac operator on the manifold with line element 
dE|,_2 [M]. From formula we obtain that the spinor 
ip of expression pTjl satisfies 



2D 



«(T3'«]V'2D ® X- 



(23) 



Thus in the Z?-dimensional spherically symmetric space- 
time with line element the Dirac equation ([T|) reduces 
to 

^2dV'2_D = ('««(T3 - imHl2)ll)2D- (24) 

Next, from the two-dimensional line element dsjD of Eq. 
(|20|) we define the line element ds^D t>y 

2 / ^2 \ / 7T, \ 2 



ds2D = ( § 



1 






' 









n- 








and usiu] 



,2 G\ 2\ F 
F2 / I iZ 



d?^ 



(25) 



(T2 


















, (73 = ( 


. 





(18) 



i'2D of the spacetime with line element ds2i5 satisfies the 
equation 

F 



that is, (Ti, (T2, and are the Pauli matrices. 

Using the basis of one-forms (|16p and the representation 
for gamma matrices given in formulae p7|) . in the space- 
time with line element ds^ of formula ([15]) we find that the 
-D-dimensional Dirac operator f/ become^ [11], [15], [16], 

m 



V2D^2D 



H 



(26) 



V = 7*Vt + Y^r + 7^Vi -I- • • • -I- 7 



Vc-2 



= (fvr^+rvr^)®i2(. 
+(T3®(yvr'^ 

l2(D-2)/2 



(19) 



where 'P 2D stands for Dirac operator on the two-dimen- 
sional spacetime with line element (is\jj defined in formula 
([251). 

Taking the variable y as 

dy ^ G 
dr F' 
we find that Eq. ([^ become^ 



(27) 



D-2y(£l-2)^ 



' 2D 



■•+7 ^ D-2 



where V^/''' and V^^ stand for covariant derivatives 
of a spinor in two and [D — 2) dimensions respectively, 7* 
and Y are a representation of the gamma matrices in two 
dimensions, 'P 2D ^s the Dirac operator on the two-dimen- 
sional spacetime whose line element is 



(7*5* + iydy)i^2D = ^(««:(T3 - imHl2)^2D. (28) 

where 7* and 7^ are a representation of the gamma matri- 
ces in two spacetime dimensions. 

Here for two-dimensional gamma matrices 7* and 7^ we 
use the representation [11], [15], [16], [17] 



7* = 



7^ = 







(29) 



ds2 



2D 



,2 G2 „ 

fiH dr 

^2 ^2 



(20) 



and ^^j^ is the Dirac operator on the {D — 2)-dimcnsional 
submanifold with line element dE^ ^ and with signature 
(+,...,+). 

For many relevant spacetimes dS]|j_2 is the line element 
of a (Z? — 2)-dimensional sphere, but this is not the only 



to find that in the Z3-dimensional spherically symmetric 
spacetime with line element ([IB]) the Dirac equation ([T]) 
reduces to 



dt^2 



G 

F , 



dr1p2 



F 

JK— — ifiF ) i/j 



F 

IK— + i^iF ] -02, 



(30) 



^In formulae ||19|I there is no sum on the repeated indices t, r, and 
D -1. 



^Notice that in Eq. I|28| l there is no sum on the repeated indices 
t and y. 
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where the functions ipi and "02 are the components of the 
two-spinor ip2D, that is 



ground is [40] 



'2D 



^2 



(31) 



Thus we get that in the ZJ-dimensional spherically sym- 
metric spacetime with line element (jl3p the Dirac equa- 
tion IT]) reduces to the pair of coupled partial differential 
equations in the variables t and r given in Eq. (j30p . This 
system of two coupled partial differential equations in two 
variables is the main result that we state in this section 
and it was previously obtained in Refs. [11], [T2] . 

It is convenient to notice that in _D-dimensional de Sitter 
spacetime, whose line element in static coordinates takes 
the form (fT3|) . when we write Eqs. ([30]) for this spacetime 
we obtain Eqs. (11) of Ref. ^33, that we previously get by 
using the results of Refs. [H] to reduce the Dirac equation. 

Although we study spherically symmetric backgrounds 
in D > 4 dimensions, we think that the results obtained in 
this section also are valid in three-dimensional spacetimes 
whose metric can take the form (1131) . We cannot compare 
in straightforward way the reduced system of partial dif- 
ferential equations presented in this section with that of 
Ref. [45] because in the previous reference a different basis 
of one-forms was chosen. 



3 Quasinormal modes of the Dirac 
field in the D-dimensional Nariai 
spacetime 

As we previously mentioned, the QNF are complex quanti- 
ties that depend on the physical parameters of a spacetime. 
Thus if we know the QNF we can infer the values of several 
physically relevant quantities of the spacetime. For many 
relevant backgrounds, for example, Schwarzschild and Kerr 
black holes, it is not possible to calculate the values of their 
QNF in exact form, we must use approximate or numerical 
methods [3]. 

Nevertheless, at present time we know many higher and 
lower dimensional spacetimes whose QNF were computed 
exactly in Refs. [21]- [39]. The systems that allow exact 
computations of some physical parameters have the ad- 
vantage that we can analyze in more detail their proper- 
ties and verify in a simple setting some predictions of the 
physical theories. Hence we believe that these examples 
deserve a detailed study. Doubtless these models will be 
useful in future research. 

As an elementary application of the coupled system of 
partial differential equations (l30|) for Dirac field propagat- 
ing in the Z)-dimensional spherically symmetric spacetimes 
that we present in the previous section, here we exactly 
compute the QNF for this field in the D-dimensional Nar- 
iai spacetime which is a simple vacuum solution of the 
Einstein equations with positive cosmological constant. 

The line element of the I?-dimensional Nariai back- 



ds^ = {l-(Tr'^)dt^ - 



dr2 



(l-(7r2) 



2\ fl^dS]|)_2i 



where 



a = {D-l)A, 



(32) 



(33) 



dE|)_2 is the line element of a {D — 2)-dimensional unit 
sphere, the constant is equal to 



(^-3) 
(i?-l)A' 



(34) 



and the constant A is related to the cosmological constant. 
If cr > then the metric ([5^ has two cosmological horizons 
at [40] 

r = ±^. (35) 
In the following we assume that the radial coordinate r e 

(-i/V^,+i/V^). 

We note that the ZJ-dimensional Nariai spacetime (|32|) 
has the following features [40]: (a) it has a geometry 
52 X S^~2, where 5*2 stands for two-dimensional de Sit- 
ter spacetime and S^~^ denotes the (D — 2)-dimensional 
sphere, (b) it is spherically symmetric, homogeneous, and 
locally static, (c) it is geodesically complete. Owing to the 
i'-dimensional Nariai spacetime p2p is spherically sym- 
metric, its metric can be written in the form (jl3p with the 
functions F, G, and H equal to 



,2x1/2 



H 



(36) 



We define the QNMs of the Nariai spacetime as the 
modes that are purely outgoing near both horizons [38], 
[39] . We also notice that the results of this section are 
an extension of those already published in the previous 
two references for coupled gravitational and electromag- 
netic perturbations, Klein-Gordon fields and tensor type 
gravitational perturbations. 

To compute the QNF of the uncharged Dirac field that 
is propagating in the Z3-dimensional Nariai spacetime ([32|), 
we first write in this spacetime the reduced system of par- 
tial differential equations for Dirac equation in D-di- 
mensional spherically symmetric spacetimes. We get the 
following system of partial differential equations 



dti^2 - (1 - ar^)dr^J2 = (1 - (Tr2)l/2 ( _ _ ) 



(37) 



dtipi + (1 - (Tr'^)^r^^Jl = -(1 -crr^)!/^ ( ^ + ) i/; 



where k are the eigenvalues of the Dirac operator on the 
ZJ-dimensional sphere, that is, k — zti{l + {D — 2)/2) and 

/ = o,i,2,..., m- 

If we take the components ipi and ^2 of the two spinor 
-0213 of formula ([5T|) as 



01 = i?i(r)e-'"*, 
^2 = i?2(r)e-'"*, 



(38) 
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then Eqs. p7|) transforms into the coupled system of ordi- 
nary differential equations 



where 



or 



2 ARi 



dr 



( 



where we define the quantity K = — k. Moreover, defining 
the following quantities z — ^/ar^ uj — aj/y/a, rh — to/^/ct, 
and A = iK/a^/a, we find that Eqs. ([39| become 

(1 _ + iuR2 = (1 - ^2)i/2(^ + irn)Ri, 

az 

(1 _ z^)^ - = (1 - z2)l/2(^ _ (40) 

dz 

Also notice that z & (—1,1). 

Next, we define (as in Chandrasekhar book's [4j) 



— arctan 



(41) 



that is 



A = VA2 +m? cos{e), 
rh — \/ + rh? sin(6'), 

and taking 

i?i = e-''/^R^, 
R2 - o^"^R2. 
we see that Eqs. pO)) reduce to 



(42) 



(43) 



(1-.^) 



dz 
2^d^i 
"d7 



= (1 - z2)1/2q,^^2, (44) 



where oat = ^/X^ + rr?. 

From the previous equations we get that the decoupled 
ordinary differential equations for the functions Ri and i?2 
are equal to 

d^^i z dRi (d)^ — iLuz)Ri aj^Ri 



dz^ 1 — z^ dz 



+ 



(1-Z2)2 



1-Z2 



d^.R2 z dR2 {Lb^ + iu}z)R2 a%R2 _ 

dz2 1 - z2 dz (1 - Z2)2 1 _ z2 " ^ ' 

To solve Eqs. (|45l) we make the change of variable 



(46) 



with y Cz (0, 1) and the ansatz 

Ri^{l-yf'y'''Si{y), 

i?2 = (l-y)^^2/^^52(y), (47) 



Bi 



^2 = 



C2 = 



2 ' 2 



' 2 



2 
2 



" 2 
' 2 



(48) 



to find that the functions Si{y) and S2{y) must be solu- 
tions of the hypergeometric differential equation 46J 

yi^-y)^ + i^-i^ + b + ^)y)^ - = O' (49) 

with parameters (the lower indices 1 or 2 determine if the 
parameter correspond to the function 5'i (y) or to the func- 
tion S2{y)) 



(50) 



ai 


= Bi4 


-Cl 


+ ia^ 


bi 


= Si4 




— iUN 


Cl 


= 2Ci 






02 


-^2 4 


-C2 


+ iapf 


b2 


= ^2 4 


-C2 




C2 


= 2C2 


+ i 





In the following we study the function i?i (we obtain 
similar results for the function i?2)- Also we take the quan- 
tities Bi and Ci as Bi = —iuj/2 and Ci = iuj/2. From the 
previous results we see that if the parameter ci is not an 
integer, then we write the function Ri as [33] 

i?i-e-^«/2(l-y)--/2 [Biy'^/\Fi{ai,bi;ci;y) (51) 
+Ei yi/2-«-/22^^(a^ - Cl + 1, 6i - Cl + 1; 2 - ci; y)] , 

where Di and Ei are constants. 

At this point we note that the tortoise coordinate for 
the D-dimensional Nariai spacetime is equal to [38], [39], 

m 



dr 



1 



:arctanh(z). 



1 — 0^2 Y 

where x S (— oo, +oo) and from expression (|46|) we get 

3^^^, and 



(52) 



-oo, 



as x^+oo, l-?/«e~2v^^. (53) 

Thus near the horizon at r = (that is as x — > 

-00) the function i?i ((5T|) behaves as 



i)ie" 



Cie"'""=e^^. 



(54) 
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In the ZJ-dimensional Nariai spacetime to satisfy the 
QNMs boundary condition near r — — that is, the 

function i?i behaves as exp {^iujx) if the tortoise coordi- 
nate goes to minus infinity x — > — oo, we take Di = in 
formula ([51)1 and therefore the function i?i becomes 

X 2i^i(ai - ci + 1, 61 - ci + 1; 2 - ci;2/) (55) 

= Eie-*«/2(i - j/)--/V/'-^''/'2J^i(«i,/5i;7i;2/)- 

We recall that if the quantity 7 — a — /3 is not an integer 
then the hypergeometric function 2Fi{a, P;j;u) satisfies 

m 



2Fi(a,/3;7;w) 



r(7)r(7 



r(7 - «)r(7 - p) 

, r(7)r(a + /3-7) 



(56) 



r(a)r(/3) 
X 2-^1(7 - "i7 - /3;7 + 1 



(1 _u)-'-"-'3 



Thus if the quantity 71 — ai — /3i is not an integer, then 
using formula (j56|) we write the radial function ()55p as 



r(7i)r(7i - cti - /3i) 

r(7i - ai)r(7i - /3i) 



X (1-2/) 



-id 12 



Fi (ai , /3i ; ai + /?! + 1 - 71 ; 1 - 2;) 



r(7i)r(ai+/3i-7i) 

r(«i)r(/3i) 

X2-Fi(7i - ai,7i - /3i;7i 



(1 - zj)l/2+^^/2 (57) 

1 - q;i - 1 - ?/)] . 



Therefore as a; ^ +00, taking into account expressions 
(j53p . we see that the function i?i is approximately equal 
to 



r(7i)r(7i - Q!i - ^1) 

r(7i - Q;i)r(7i - ^1) 
r(7i)r(ai + /3i - 71) 



r(«i)r(/3i) 



(58) 



The boundary condition for QNMs of the £)-dimension- 
al Nariai spacetime imposes that the function i?i behaves 
in the form exp (jwx) as a; ^ +00. Thus to satisfy the 
boundary condition of the QNMs for Nariai spacetime near 
the horizon at r = +1/ ^Ja we must cancel the second term 
in formula (|58p . One way is to exploit the zeros of the 
terms l/r(x) which are located at a; = — n, = 0, 1, 2, . . . . 
Hence to satisfy the boundary condition near the horizon 
at r = +\/ ^/a we must impose the condition 



ai = — rt. 



/3i = -71, 



(59) 



which imply that the QNF of the Dirac field in ZJ-dinien- 
sional Nariai spacetime are determined by the expression 



1 

UJ — iofAT — i \ n + — 



(60) 



A similar computation for radial function R2 also yields 
the QNF of formula pi|. 



Taking into account that 



^,^(2l + D-2l^ (61) 



we find that in the D-dimensional Nariai spacetime ((32l 
the QNF of the Dirac field are equal to 



UJ 



{21 + D - 2) 
4a2cr 



iJ(T n + 



1 



(62) 



In our notation the previously calculated QNF for Klein- 
Gordon field and tensor type gravitational perturbation 
are written as [39| 



[lil + D-?,)-l^f' 



iJa n + 



(63) 



We point out that for / = the QNF ^ of the Klein- 
Gordon field and tensor type gravitational perturbation 
are purely imaginary. This fact was not noted in Ref . )39j . 

When in the _D-dimensional charged Nariai spacetime 
studied in Ref. [55] we take the electric charge of the space- 
time equal to zero, for electromagnetic and gravitational 
perturbations of vector type we get their QNF are equal 
to 



(/(? + i?-3) + 2i2_i3)i/2 1 



UJy=± 



(64) 



whereas for electromagnetic and gravitational perturba- 
tions of scalar type their QNF are 



= ± 



1/2 



iy/a \ n + 



1 



(?(/ + D-3)-^^)i/2 



— i\/(T n 



(65) 



Notice that in Refs. [351, [SH], to calculate the QNF ([55)) . 
([M|) . and ([^5]) of the D-dimensional Nariai spacetime the 
result for QNF of Poschl- Teller potential was used. The 
QNF of this potential were previously computed in the 
paper by Ferrari and Mashhoon [47] , 

It is convenient to observe that for identical values of 
D, n, and A the imaginary part of the QNF ([62]) for Dirac 
field is identical to the imaginary part of the QNF for grav- 
itational and electromagnetic perturbations of vector type 
and scalar type (formulae ([M|) and ([SS]) ). and for Klein- 
Gordon fields and tensor type gravitational perturbations 
(formula (|63p ) already computed in Refs. [38], [39]. Thus 
the decay time r — l/Im(cL;) is the same for fermion and 
boson fields studied here and in Refs. [38], (3^. Also we 
point out that for QNMs of these fields the decay time 
does not depend on the angular momentum number I. 
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For Dirac field the real part of the QNF ((62)) show some 
differences with respect to real part of QNF for Klein- 
Gordon, electromagnetic, and gravitational perturbations 
dMl), (IMl), and dnni). We think that the source of these dif- 
ferences is that in the present section we study the Dirac 
field, whereas in Refs. [38] , [39] the fields studied are mass- 
less boson fields. 

From formula ([S^ . we obtain that for massless Dirac 
field (Weyl field) its QNF are equal to 



21 + D-2 
2a 



(66) 



Even in this case the real part of the QNF (|66l) shows a 
different dependence on the angular momentum number 
that the real part of the QNF for Klein-Gordon, electro- 
magnetic, and gravitational perturbations. The imaginary 
part is equal to that of the Dirac field. Thus for Dirac field 
in the Z)-dimensional Nariai spacetime the decay times of 
its QNMs do not depend on the mass. 

Taking into account the results for QNF (|^ . (IM)) . 
(|65|) . and (|66|) of the £)-dimensional Nariai spacetime, we 
note that their real and imaginary parts show an explicit 
dependence on the dimension D of the spacetime (the 
imaginary part through the parameter a of formula ([33])). 
For all the fields whose QNF have been calculated and 
for a given mode number n we infer that the decay time 
decreases as the dimension of the spacetime D increases. 
The dependence of the oscillation frequency on the space- 
time dimension D is more complicated, it first decreases 
and then increases as D increases. Moreover for / > 3 and 
for the same values of -D, I, and n the oscillation frequen- 
cies of the massless Dirac field (l66t are greater than the 
oscillation frequencies for Klein-Gordon, electromagnetic, 
and gravitational perturbations of formulae ([63]) . ([64]) . and 
([65]). 

As the harmonic time dependence is of the form 
exp(— iwi) (see formulae ([55]) '). in order to have stable 
QNMs we need that Im(cj) < 0. We notice that for QNF of 
formula Im(tj) < 0, thus the QNMs of the Dirac field 
decay in time. Also for Klein-Gordon, electromagnetic, 
and gravitational perturbations a similar result is true (see 
formulae (I63|) . ([64]) . and ([65]) ). Therefore we assert that 
the D-dimcnsional Nariai spacetime is perturbatively sta- 
ble under the propagation of these classical fields. 

To finish this section, we note that in Refs. [38], [39] 
is shown that the radial differential equations for Klein- 
Gordon fields, tensor type gravitational perturbations, and 
coupled electromagnetic and gravitational perturbations 
reduce to Schrodinger type equations with a Poschl- Teller 
potential of the form 



U 



cosh {-sf^x) 



(67) 



where the value of the constant U depends on the pertur- 
bation type (see Refs. [38], [39]). 

Following the procedure of Chapter 10 in Ref. [1] we 
transform Eqs. ([55]) into a pair of Schrodinger type equa- 



tions with potentials equal to 



V±{x) 



cosh [y/ax) 



(68) 



± 



/ (21+D-2) 



1/2 



a sinh(Y^2;) 



cosh {^/ax) 



We observe that these potentials are of Morse type (see 
Table I of Ref. [IS]). 

For plots of the potentials (|67p and ([68| see Figures 1 
and 2. We note that for identical values of D, A, and I the 
shape of the potentials (|67p and ([55]) is similar, only ob- 
serve that the height of Poschl- Teller potential is smaller 
than the height of Morse potentials. Notice that in Figure 
1 is plotted the Poschl- Teller potential corresponding to 
the QNF ujy, but for Poschl- Teller potentials ([67]) corre- 
sponding to the QNF ujy, ujg, w^, and ujkg also is true 
that for identical values of D, A, and I their height is 
smaller than the height of Morse potential ([S5]) . 



T 




P 
/ 20 




/ 




/ 




/ 5 





Figure 1: Plot of Poschl- Teller potential V of Eq. ([67]) 
corresponding to cuy, where we take D = 5, A = 1, and 
1 = 3. 




Figure 2: Plot of Morse potentials V+ (solid line) and 
(dashed fine) of Eq. ([68]). where we take m = 0, D = 5, 
A = 1, and / = 3. 



4 Concluding remarks 

The _D-dimcnsional Nariai spacetime studied in Section[3]is 
uncharged. Notice that there is a charged generalization 
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of the D-dimensional Nariai spacetime, only we need to 
replace the quantity a of formula ((33l) by 0, [40] 
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aQ = {D-l)A- 



(D - 3)2g2 



a 



2{D-2) 
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charged case by ag which is a solution to the equation 



{D-3) 
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2(D-2) 
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where Q is related to the electric charge of the spacetime. 

Supported in our mathematical analysis of the problem 
for the uncharged Nariai spacetime, we assert that in the 
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